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ABSTRACT 



Context. A 321.5 s modulation was discovered in 1999 in the X-ray light curve of HM Cnc. In 2001 and 2002, optical photometric 
and spectroscopic observations revealed that HM Cnc is a very blue object with no intrinsic absorptions but broad (FWHM~1500 km 
^rfH ' s"') low equivalent width emission lines (EW~ 1 -r 6A), which were first identified with the Hell Pickering series. The combination 

I of X-ray and optical observations pictures HM Cnc as a double degenerate binary hosting two white dwarfs, and possibly being the 

\^ , shortest orbital period binary discovered so far. 

Aims. The present work is aimed at studying the orbital motion of the two components by following the variations of the shape, 
centroid and intensity of the emission lines through the orbit. 

Methods. In February 2007, we carried out the first phase resolved optical spectroscopic study with the VLT/F0RS2 in the High Time 
^ , Resolution (HIT) mode, yielding five phase bins in the 321 s modulation. 

\^ • Results. Despite the low SNR, the data show that the intensity of the three most prominent emission lines, already detected in 2001, 

' varies with the phase. These lines are detected at phases 0.2- 0.6 where the optical emission peaks, and marginally detected or not 

I detected at all elsewhere. Moreover, the FWHM of the emission lines in the phase resolved spectra is smaller, by almost a factor 2, 

than that in the the phase-averaged 2001 spectrum. 

Conclusions. Our results are consistent with both the pulsed optical component and emission lines originating in the same region 
which we identify with the irradiated surface of the secondary. Moreover, we note that the EWs of the emission lines might be up to 
^ ' -15 -r -25 A, larger than thought before; these values are more similar to those detected in cataclysmic variables. All the findings 

^ : further confirm that the 321s modulation observed in HM Cnc is the orbital period of the system, the shortest known to date. 

00 . 
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(3 1 . Introduction In 1999, Israel et al. (1999, hereafter 199) discovered an X- 

I ray emitting source with a strong 321.5 s. Deep TNG and VLT 

observations were carried out to study the optical counterpart 

_ , , , , . . .T^r^x, N , , • of HM Cnc and its optical spectral features. HM Cnc is a blue 

Double degenerate binaries (DDEs) are close bmary star systems oi i /d » u • icw i » »u 

. . ° J. , . , ^ /„,T^s i,,,^ , r. V=21.1 (B=20.7) star showing a ~15% pulsed fraction at the 

consisting of two white dwarfs (WD) or a WD and a He star. ^ ^„ • j /t i ^ i ^^vrv^ u rm i 

° , . . ^ . ,. , ~32 1. 5 s X-ray period (Israel et al. 2002, hereafter 102; see also 

,^ . DDBs are expected to originate from main sequence stars which „ , , Inn^ ^ rr, i j v i j- i 

5-H , . ^ , ° , , ^ Ramsay et al. 2002a). The optical and X-ray pulsations displays 

' have experienced two common envelope phases, eventually ex- u cp ^ c n on n oc /r> » i onrn\ • j- . 

... . , ,, J, , • • , • a phase offset of ~0. 20-0. 25 (Barros et al. 2007) indicating that 

posing the CO or He cores of the original stars (for a review see • • . ^ . jvcr . • t-u i j i 

i,, ° . ^, , • , . , they originate from two different regions. The optical modula- 

Warner 1995). Among DDBs are AM CVns, systems in which , • . • . -.u u-i i j f 

' ,f,_ , , ■' , tion also varied in shape with time, while its pulsed fraction 

either a low mass He WD or a low mass non degenerate He burn- , , , , j wu • i 

^„ . „ , , , , ^ ° , . was clearly wavelength dependent (being larger at longer wave- 

ing star fills its Roche lobe and transfers matter to the companion , \ -ri, . f i ii\ ^ /-^ • u . • j u ui 

f , , N , ,. • , • , , lengths). The spectrum of HM Cnc IS characterized by a blue 

(Nelemans et al. 2001a). AM CVn binaries are characterized by ~. .J • • u i- /Tm\ u . u j 

^ , . , ... , ,„„ ^^^^ , • , , continuum with no intrinsic absorption lines (102), but broad 

orbital periods m the 600-3000 s range and, in the case of non de- , , _is , ■ , •,, ,t^„, , ^ 

generate donors, they evolve from long to short orbital periods (FWHM~1500 km s '), low equivalent width (EW -1-6 A 

Ti n u »u • • 1 • J f in • N • • emission lines, which Were first identified as Imes from the Hell 

(till they reach the minimum orbital period of ~ 10 mm), simi- „. , . • . , n „ , • . 

, , . . , • • , , .^x, s rj,, , . ■ , Pickering series. All these findings are fully compatible with 

larly to cataclysmic variables (CVs). They evolve toward longer . . ?, . , r^rir, , 

orbital periods, otherwise (Nelemans et al. 2001a and reference C"'^ being an X-ray emitting. He-rich DDB, possibly a pro- 

.■ • c ^ .. A ■ » e genitor of the AM CVn systems, with the shortest orbital pe- 

therein). Gravitational radiation waves drive the mass transfer ° , , , ..^ . „,,, ^ , 

■ ■ , , , ■ .... ... . , . f . nod ever recorded (102). Due to the ditterence in EWs measured 

and double-peaked emission lines in the optical spectra of most , , , , , , , , tt . • 

\\/< r<\r . . £ . f t- A- \ ^ A- t- between the odd- and even-terms of the Hell Pickering series, 

AM CVns testifies to the presence of an accretion disk mediating . , , ^ , , , , 

a e e » » »u »u It has been proposed that HM Cnc may also be hydrogen-nch 

the flow of matter from one star to the other. ^ ^ ^ i s= 



X 
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(Steiper et al. 2005 and references therein). Should this be con- 
firmed, it would be diffcult to relate HM Cnc with AM CVn 
systems which show only He emission/absorption lines in their 
spectra. Thanks to the 200 1-2006 photometric monitoring, Israel 
et al. (2004) found a very accurate P and P phase coherent solu- 
tion, with P - -1.1x10"^ s/yr implying that the orbit is decay- 
ing. The optical solution was sufficiently accurate that it could 
be extended backward in time so as to find a unique phase- 
connected solution encompassing ROSAT 1994 and 1995 obser- 
vations. This yielded a 10-years baseline P-P coherent solution 
giving P=321.53038(2) s and P = -3.661(5)xlO-" s s"' (for 
more details see Israel et al. 2004; see also Strohmayer 2005 and 
Barros et al. 2007). Yet, the nature of the X-ray emission de- 
tected from HM Cnc and its twin source V407 Vul is still under 
debate. A number of models have been proposed (for a review 
see Cropper et al. 2004). Among these, is the DDB model with 
mass transfer proposed in two flavors: with a magnetic primary 
(polar-like. Cropper et al. 1998) and with a non-magnetic accre- 
tor (Algol-like, Mai'sh & Steeghs 2002, Ramsay et al. 2002b). In 
the latter model, the reason why a disk does not form is because 
the minimum distance of the gas stream from the center of mass 
of the system is smaller than the size of the accretor, resulting in 
the stream hitting the surface of the accretor The unipolar induc- 
tor model, UIM, is a second accredited model and, according to 
it, the secondary star does not fill its Roche-lobe but, by moving 
across the primary star magnetic field, induces an electric current 
which heats up the B-field line footprints on the primary white 
dwarf. No mass transfer is expected within the UIM. A similar 
mechanism is likely responsible for the discovered Jupiter-Io in- 
teraction (Wu et al. 2002; Clarke et al. 1996; Dall'Osso et al. 
2006, 2007). In the DDB scenario the X-ray emitting region is 
thought to be "self-eclipsed" by the primary object giving rise 
to the sharp and large amplitude X-ray modulation. In all mod- 
els the source is expected to be one of the strongest sources of 
low frequency gravitational wave radiation, easily detectable in 
the future by the LISA mission (Nelemans et al. 2001b and ref- 
erences therein). Within the DDB scenario, the negative value 
of the period derivative suggests that matter transfer is not vi- 
able, therefore disfavoring the Polar-like and Algol-like scenar- 
ios, and the fact that the GW emission is driving the orbital pe- 
riod decay. Moreover, the P value is slightly large for the grav- 
itational wave emission mechanism as the only channel govern- 
ing the period evolution and we cannot currently exclude that 
magnetic stresses are also present in the system (Dall'Osso et 
al. 2006, 2007). However, Deloye and Taam (2006) have revised 
AM CVn formation and evolution, showing that fully degener- 
ate donors and the accretion scenario can easily produce systems 
having the orbital periods of HM Cnc and V407 Vul at early con- 
tact. In addition, they also revise the duration of the phase char- 
acterized by mass transfer and negative P, determining values of 
Tj^ which are comparable with those predicted by the UIM. An 
alternative way around the negative P during accretion has been 
proposed by D' Antona et al. (2006), who suggest that HM Cnc is 
an accreting DDB, which is characterized by p-p hydrogen shell 
burning in the donor, seen during a mass transfer phase. The nu- 
clear reactions on the shell cause the secondary star to shrink 
in response to the mass loss. D'Antona et al. (2006) obtain a 
period derivative which is correct in both sign and magnitude. 
This scenario also accounts for the possible presence of H in 
the spectrum of HM Cnc. The lack of (relatively large) optical 
polarization argues against the Polar-like scenario. We will not 
consider this model further in this paper (see Reinsch et al. 2004; 
Israel et al. 2004). Recently, Ramsay et al. (2007) reported the 
detection of transient radio emission from HM Cnc, which can 



be accounted for only with maser emission in the framework of 
the UIM model. 

In 2006 we requested and obtained VLT-FORS2 observa- 
tions in a new observing mode, aimed at carrying out the first 
ever phase-resolved spectroscopic study of the 5.4 min modula- 
tion of HM Cnc. Our main goals included the disentangling of 
the different optical components (pulsed vs un-pulsed), the iden- 
tification of the emission line forming region(s) and the detection 
of the emission lines Doppler modulation to provide the ultimate 
proof of HM Cnc binary nature with orbital period OP=5.4 min. 
Here, we report on the results we obtained from this new spec- 
troscopic campaign. In section 2 we describe the observation and 
the FORS2 HIT-S mode. In section 3 the obtained results, while 
in section 4 we report our conclusions. 

2. Observation and data reduction 

The data were collected during four half nights at UT1+FORS2 
in High Time Resolution spectroscopy (HIT-S) mod^l This ob- 
serving mode allows fast time resolved spectroscopy with little 
overheads by shifting the charges on the detector while expos- 
ing. Using a mask on which a small slit has been cut, only a 
small part of the detector get actually exposed to light. The ex- 
posed part is then shifted under the masked area while an new 
spectrum is taken. The sequence continues till the whole CCD 
has been used-up and it has to be read-out. By using a mask with 
slit 1.0" wide and 10" long, in 4 nights, we collected 15 frames 
containing 41 spectra each. Each spectrum corresponds to 64 s 
of integration, hence the total exposure time per frame adds-up 
to 2624 s (-44 min). 

The adopted instrument setup was grism 600B covering the 
3550-6600A wavelength range. A gap of ~10A at about 5000A 
was present due to the particular orientation of the grism in the 
HIT mod^ The resulting spectral resolution was ~5A. 

The data produced within the FORS HIT-S mode looks sim- 
ilar to multi object spectroscopy data, being they a collection of 
spectra in a single frame. However, as we experienced during 
the data reduction process, the handling of HIT-S mode frames 
requires particular attention. In particular, we have verified that 
the standard frame pre-processing (bias subtraction, flat correc- 
tion etc) is not correct as it introduces noise. The HIT-S mode 
frames present a non uniform and variable bias level which in- 
creases in the direction of the charge shift. In the majority of 
our science frames we observed a sloped bias with an increase 
of ~5-7 ADUs, or 2-4%, from left to right. In three frames the 
"bias slope" was steeper with the bias level increasing by 14% 
or up to 250 ADUs across the detector We discarded those three 
frames in order to work with an "homogeneous" data set and re- 
duced the remaining 'good' frames subtracting, in place of the 
bias, a third order polynom fitted along the unexposed rows of 
the detector 

We also discarded the flat fields for the following two rea- 
sons. First, their showed yet a different and variable "bias level". 
Second and most important, the flat field and our science data 
are characterized by completely different counts level regimes 
so that the first does not correctly remove the pixel-to-pixel vari- 
ations of the second. Due to the relatively high spectral reso- 
lution, the faintness of our target and the short exposure time, 
the signal of the object-Hsky was only a few ADUs above the 

' http://www.eso.org/instruments/forsl/ for more information. 

' The orientation is 90° with respect to the standard FORS configura- 
tion. Hence, the spectra are no longer parallel to the CCD-mosaic rows 
but to their columns. 
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Fig. 1. The five 2007 VLT-FORS2 phase-resolved spectra of HM Cnc are shown plus, at the top, the 2001 FORS 1 and 2007 FORS2 
phase-averaged spectrum (adapted from 102). All the spectra are normalized to one and arbitrarily shifted for clarity. The numbers 
on the left of the spectra mark the center of the 0.2-long phase interval. Vertical dashed lines mark the Hell 4100.0 A, Hell 4338.7 A, 
and Hell 4859.3 A. The spectra were re-binned at the spectral resolution imposed by the slit width (4.5 A). 



bias level. In low count level regimes of 10-50 e /pix, CCDs 
suffer from deferred charges or charge traps (Mc Lean 1997). 
Hence, charge traps are most visible in the case of short ex- 
posures or spectroscopic observations (for which the dispersed 
background is below lOOe ). Our dispersed signal is of the order 
of lOe , only. Delayed release of trapped charges can occur in 
time ranges from a few second to hours. The net result will be 



1-2 pixels wide dark (or bright) short strips in the direction of 
the charges tansfer and an overall increase of the noise. We note 
that in our case these features are of the same order of the noise 
and therefore are not visible in the single spectra but only in the 
final median/average spectra for which the noise is significantly 
reduced. We have verified that this is the case by reducing the 
same data sets following two different strategies (see below) and 
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comparing the results. We have identified at least 6 CCD rows 
afi'ected by charge-traps. The release of the trapped charges oc- 
curs randomly and cannot produce spurious emission features as 
would not always appear in the same position on the detector. 
Hence, they are removed when the median of the phase-bins is 
computed. 

We continued the data reduction following two different 
strategies both as a cross check and to pinpoint the best data re- 
duction procedure. In both cases we phased each single spectrum 
of each frame and defined 5 phase bins centered at phase 0.1,0.3, 
0.5, 0.7 and 0.9 (see below for the absolute time and phase as- 
signment procedure). In one case we first sliced each frame in 
41 subframes. The 2D slices were then stacked and their median 
computed in one of the five phase-bins which where eventually 
extracted. All the spectra within a given shce were shifted to a 
common value in the spatial direction, in order to avoid smearing 
of the signal because of the limited centering accuracy during the 
acquisition. The 5 median phase-bin spectra were wavelength 
calibrated in 2D before extraction. 

In the second case all the spectra were extracted from the pre- 
processed 2D frame, after having build a reference trace function 
on the median spectrum of the given night. The single extracted 
spectra were wavelength calibrated and combined via median 
in the 5 bins. In both cases we did not combined all the pre- 
processed frames but performed a visual inspection discarding 
those spectra that had a particularly bad SNR because of the poor 
sky conditions. As a result of this selection 75, 79, 75, 73 and 69 
spectra were combined in bin 1,2,3,4, and 5, respectively, out of 
the original ~ 123 spectra per phase bin. We verified that this pro- 
cedure is equivalent , if not better, to combining all the spectra 
(good and bad ones) with different weights (e.g. accounting for 
the slit losses due to variable seeing conditions). 

2.1. Timing and pliasing each spectrum 

Particular attention was payed to the time shift introduced by the 
(solar system) barycentric correction which is of the order of few 
seconds within each night and up to ~30 s (~0.1 in phase) over 
the whole observing run. We computed, a priori, the barycentric 
correction to be applied to each frame starting times. This dras- 
tically reduced the overlapping phase interval among phase bins 
(the overlapping being of order 50% for the uncorrected frames). 
The time accuracy of the five phase bins is of order 20-25% or 
~14 s out of the 64 s exposure time of the spectra. The ~14 s 
uncertainty is mainly due to the unfortunate use of 64 s instead 
of 64.3 s (= 321.53 s/5) as effective exposure time for spectra. 
The ephemeris used for the phase assignment are those of the 
coherent timing solution by Israel et al. (2004), and extrapolated 
to MJD 54144.0 (the first observing night of the run). The period 
extrapolated to MJD 54144.0 is P=321.52555(3) s. Note that we 
did not correct for the P variation during the observing run base- 
line, since it accounts for up to a maximum of 2x10"^ s within 
five nights, well below the 14 s phase accuracy. 



3. Phase-resolved spectroscopy 

Given the relatively poor statistics, we did not attempt to study 
the continuum of each phase-resolved spectrum, nor we looked 
for new and/or unexpected features (whether in emission or 
in absorptions). We rather focused on the normalized phase- 
resolved bin-spectra searching for the strongest emission lines 
already detected in 2001 (see Figure 1). A visual inspection re- 
sulted in the possible identification of three emission lines at 



about 4100 A, 4340 A, and 4860 A in the spectrum centered 
at phase 0.5. Therefore we fitted the normalized spectrum with 
three Gaussians centered at the three wavelengths reported above 
and measured for each of them the intensity, FWHM, Ac, and 
EW. In the phase 0.5 spectrum the detection of these lines is con- 
firmed at a ^3cr level and the best fit parameter are reported in 
Table 1. By applying the same analysis to the remaining spectra, 
we obtained only marginal detection (~2-3cr level) of the same 
lines in the spectrum centered at phase 0.3 and no-detection in 
the spectra centered at phase 0.7 and 0.9. In the phase 0.1 spec- 
trum we detect only the Hell /14859 at ~3cr level (see Table I). 
At last, we also fitted the same three lines both in the phase- 
averaged spectum obtained by combining the 5 phase-bin spec- 
tra, and in the VLT-FORSl 2001 spectrum (see 102 for details). 
While in the 2001 spectrum the detection of the emission lines 
at 4100, 4338 and 4859 A is highly significant (>6cr), in the case 
of the average 2007 spectrum the same lines are only marginally 
detected as expected when combining a small signal with no sig- 
nal. However, in the average 2007 spectrum we detect at a ~3cr 
confidence level an emission line centered at ~4700 A. This was 
identified with a blend of CIH, NIH and Hell in the 2001 spec- 
trum (Israel et al. 2002). 

Table 1 reports the best fit parameters of the detected emis- 
sion lines together with the upper limits on the intensity of the 
phase-resolved and phase-averaged spectra. The upper limits on 
the intensity of emission lines were obtained in two different 
ways. We first evaluated the significance of any additional com- 
ponent (each of the three Gaussian) to the fit of the normaUzed 
continuum (set to one by definition) applying the F-test. The F- 
test probabilities are reported in Table 1 in cr units. However, we 
note that Protassov et al. (2002) pointed out that the F-test may 
be inappropriate in some circumstances and lead to incorrect es- 
timates of the significance of a feature. In order to further in- 
vestigate this issue, we ran a Monte Carlo simulation of 100000 
spectra with only the continuum model, and the same spectral 
coverage and resolution of the spectra we have analyzed. The 
results of this simulation gave very similar values to those in- 
ferred with the F-test, confirming the confidence levels reported 
in Table 1 . Note that the above reported confidence levels refers 
to a search over the whole spectrum (the number of trial is the 
number of independent bins of wavelength); however our emis- 
sion line search has been restricted only to the strongest emission 
lines detected in 2001 with a number of trial equal to few for 
each search. Correspondingly, the the reported confidence levels 
are lower limits. Moreover, we can confidentely exclude that the 
observed line variability is originated by SNR variations (among 
phase resolved sepctra) due to the intrinsic flux variations (of the 
order of 10-15% at maximum) of the source, since the spectrum 
where we detected most of the lines (see Figure 1) is at 0.5 in 
phase, close to the optical minimum (see Figure 2). 

It is apparent that three emission lines are significantly de- 
tected at phase 0.5, marginally detected at phase 0.3 and unde- 
tected in the three remaining phase-resolved spectra (with the 
exception of the line around Hel /i4859 at phase 0.1). The com- 
puted upper limits are such that we can confidently conclude that 
the emission fines in the spectra centered a phase 0.1,0.7 and 0.9 
are either absent or considerably weaker than in the spectra cen- 
tered at phase 0.3 and 0.5. 



4. Discussion 

The VLT-F0RS2 campaign we carried out in February 2007, 
though affected by poor sky conditions (seeing in the range 1.2- 
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Fig. 2. The 1994-2002 phase coherently connected X-ray folded light curves (filled circles; 100% pulsed fraction) of HM Cnc, 
together with the VLT-TNG 2001-2004 phase connected folded optical light curves (blue asterisks). For both light curves we have 
applied the timing solution reported in section 2, after having checked that the arrival times were barycentrized and registered with 
respect to the same reference system (Terrestrial Time; see also Barros et al. 2007 for a coherent timing solution which includes the 
covariance terms of the fit ). Two orbital cycles are reported for clarity. The confidence level (in units of cr) of the three detected 
emission lines in the phase resolved VLT-FORS2 spectra are also shown (filled triangles, squares, and empty stars for 4100.0 A, 
4338.7 A, and 4859.3 A emission lines, respectively): it is evident the correlation with the optical pulsed component. 



2", and thin-thick clouds), allowed us to obtain the first phase- 
resolved spectra of HM Cnc and the following two major find- 
ings: 

- Line variability The three most intense emission lines iden- 
tified in the 2001 phase-averaged spectrum are significantly 
detected only in spectrum centered at phase 0.5, marginally 
detected in the spectrum at phase 0.3 and absent or consider- 
ably weaker in the spectra centered at phase 0. 1, 0.7 and 0.9. 
This implies that the intensity of the emission lines varies 
across the period. 

- Line width We measure narrower FWHMs for the emission 
lines (particularly in the case of lines Hell lines 4100 and 
4338A) in the phase resolved spectrum centered at phase 0.5 
than in the phase-averaged spectra taken in 2001 and 2007. 
The most straightforward interpretation is that the emission 
lines of the phase averaged spectra are Doppler broadened. 



However, the accuracy of our measurements is not sufficient 
to be conclusive about this issue. 



Additional observational constraints can be derived by com- 
paring the optical and X-ray light curves with our phase resolved 
spectra. In figure 2 we show the X-ray (black filled circles) and 
the optical (blue crosses) light curves together with the the sta- 
tistical significance of the detected emission lines (red symbols, 
different symbols are for different emission lines), as reported 
in Table 1. It is evident that the intensity of the emission lines 
peaks at the time of the X-ray maximum/pulse. However, the 
lines are likely present (absent) also at orbital phases when the 
X-ray pulse is off (on). At the same time, the emission lines re- 
gion (ELR from now on) and the source of the optical pulse are 
both off in the phase interval ~0.6-0.9, though the ELR stays 
undetected possibly till phase 1.1. 
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Table 1. Measured parameters of the emission lines detected in the phase resolved VLT-FORS2 spectra of HM Cnc. 



Line ID 




Phase resolved 






Phase averaged 


A 


0.0-0.2 


0.2-0.4 


0.4-0.6 


0.6-0.8 


0.8-1.0 


2007 


2001 


Hell 


I<0.07 


1=0 07+" 


I=0.16±0.04 


I<0.08 


I<0.07 


I=0.05±0.03 


1=0.06+0.01 


4100 




FWHM=24±7 


FWHM=16±3 






FWHM=17±f 


FWHM=22±4 






i,=4087+6 


^<.=4095±3 






i,=4097±^ 


^=4104±3 






EW=4.2±3* 


EW=6.4±2.0 






EW=2.1±f° 


EW=3. 3+0.8 






[2.80-] 


[3.90-] 






[1. 60-] 


[1.50-] 


Hell 


I<0.07 


1=0.11+0.03 


i=o.i8±^:^j 


I<0.06 


I<0.07 


I=0.09±0.03 


1=0.09+0.01 


4338.7 




FWHM=26±9 


FWHM=13±3 






FWHM=29±8 


FWHM=25±4 






^=4348+4 


^<.=4339±2 






i,=4336±3 


^r=4341±2 






EW=7.2±3.2 


EW=5.9±^J 






EW=6.5±2.8 


EW=5.6±1.1 






[3.10-] 


[3. 80-] 






[2.50-] 


[9.O0-] 


Hell 


RO.13+0.04 


1=0 06+" 


I=0.13±0.05 


I<0.05 


I<0.06 


I=0.05±0.02 


1=0.08+0.01 


4859.3 


FWHM=26±9 


FWHM=33±10 


FWHM=23±J« 






FWHM=32±12 


FWHM=39±5 




^c=4864±4 


^,=4871 ±8 


FWHM=4856±3 






/l,=4862±6 


^,=4860+2 




EW=8.5±3.9 


EW=5.0±2.2 


EW=7.5±^° 






EW=4.0±2.2 


EW=7.8+1.4 




[3.20-] 


[2.10-] 


[2.90-] 






[2.0£r] 


[8.O0-] 


blend~4700 


I<0.07 


I<0.07 


I<0.07 


I<0.07 


I<0.07 


I=0.07±0.02 


1=0.07+0.01 



FWHM=15+4 FWHM=18.4+2.5 

/l,=4696+5 ^,.=4692.5+2.5 
EW=2.6+1.0] EW=3.2+0.6 
[3.50-] [Toi 

Note - Emission line intensity (I) are in ADUs; FullWidth Half Maximum (FWHM), central wavelengths (Ac), and equivalent widths (EW) are in units of A. Intensity upper limits to 
emission lines are at lo" contidence level. The emission line detection level is reported in square bracket in units of cr (see text for more details). All the reported uncertainties are at Itr 

level. 



X-X-ray pulse, 0-optical pulse. L-emission line region 



case-A 



case-B 



case-C 




Fig. 3. A sketch of the three configurations described and dis- 
cussed in section 4. The actual position of the labels X,0 and 
L (described in the figure) on the two star components is just 
indicative. See text for mode details. 

We can imagine three different scenarios depending on the 
geometrical location of the X-ray and optical pulse and the ELR 
(see Fig[3]for a sketch): 
_ case A: 

a hot spot on the primary white dwarf is the source of the X- 
ray pulse and of the optical emission lines, while the secondary 
star, irradiated by the primary, is the source of the optical pulse 
(FigOl top panel). This is the most natural configuration if one 



considers the phase match between the X-ray maximum and the 
maximum of the emission line intensity as a starting point. The 
secondary star can then be irradiated either by the X-ray emit- 
ting hot-spot or by the hot primary white dwarf. In the former 
case, the X-ray emitting hot-spot has to be at an angle >90'' with 
the line connecting the two stars, in order to be the irradiating 
source. However, this geometry makes it difficult to account for 
the phase lag which we observe between the optical and the X- 
ray pulse (the phase lag should approach 0.5 in phase and not 
0.2-0.25, but see also Barros et al. 2007). In the case the irradi- 
ating source is the primary WD, on the other hand, there would 
be no additional constraints on the position of the X-ray hot spot 
(a part that imposed by the measured phase offset between the 
X-ray and the optical light curve). Note that, in this latter case, 
the WD has to be the source of the optical un-pulsed compo- 
nent at temperature is T>40000 K (Israel et al. 2002), in order 
to be the irradiating source. In this case-A configuration, inde- 
pendently on the irradiating source, the large observed FWHMs 
of the emission lines can be explained either by Doppler broad- 
ening by the rapidly spinning primary, close to the white dwarf 
break-up velocity (which represents a severe constraint), or by 
strong turbulence and large thermal velocities of the gas in the 
hot-spot region. The observed FWHM can be still accounted for 
by orbital Doppler smearing if the mass ratio is larger than about 
1/2. In addition, the early appearance of the ELR with respect 
to the X-ray pulse poses further constraints without an immedi- 
ate/obvious explanations. Namely, while temperature gradients 
in the hot-spot region can be invoked in order to account for 
this observation, the lower temperature region is located "down- 
stream" or beyond the X-ray emitting (impact) region; contrarily 
to what typically observed in the case of accretion induced hot- 
spots (e.g. WZ Sge in Mason et al 2000). However, an impact 
region with a down stream temperature gradient could be pro- 
duced by a fast rotating WD as explained by Marsh and Steeghs 
(2002) in the direct impact model 

_ case B: 

As in case A the X-ray hot spot is on the primary star However, 
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now, both the pulsed optical component and the ELR are on the 
secondary star (Figl3] middle panel). This configuration is well 
supported by the measured FWHM of the emission lines, which 
is consistent with Doppler broadening by the secondary star or- 
bital motion, and by the good agreement between the optical 
pulse and the visibility interval of the ELR. In this configuration, 
the phase range 0.6-0.9 match the the time of the secondary in- 
ferior conjunction when the irradiated face of the secondary star 
is not visible to the observer because of the star orbital motion 
around the binary system center of mass. Similarly to case A, the 
phase offset between the X-ray and the optical (plus ELR) pulse 
is better accounted for assuming that irradiation of the secondary 
star is due to the hot primary and not to the X-ray emitting hot- 
spot. However, within this scenario there is not an obvious ex- 
planation for the fact that the ELR and the X-ray pulse reach 
maximum at the same phase. 

We also note that the un-pulsed optical continuum accounts 
for up to 85% of the total observed B band flux, regardless of its 
origin (primary or secondary star). Therefore, in this cofiguration 
where the ELR and the optical pulsed component match (they 
both originate from the irradiated surface of the secondary), the 
equivalent width of the emission lines, after subtraction of the 
un-pulsed continuum, would be 10-25 A. These values are typi- 
cal of AM CVn systems and, more in general, cataclysmic vari- 
ables. 
_ case C: 

In this latter configuration the X-ray hot-spot, the pulsed opti- 
cal component and the ELR are all located on the primary star 
(FigE bottom panel). In order to explain the fact that the optical 
pulse and the ELR appear before the X-ray pulse we can imag- 
ine an impact region with a thermal gradient down stream due 
to the rapidly rotating underlying white dwarf (as explained by 
Marsh and Steeghs 2002 in their direct impact accretor model). 
Alternatively the three components, though all on the primary, 
might be unrelated and have different origin: e.g. the X-ray pulse 
could come from the hot-spot, while the optical pulse plus the 
ELR would come from the primary star side which is irradi- 
ated by the hot secondary, in the hypothesis that the latter is the 
source of the T>40000 K un-pulsed optical component. As in 
case A, the large FWHM could be interpreted either with tur- 
bulence and high impact or thermal velocities or with a rapidly 
spinning white dwarf primary. Similarly to case A, the observed 
FWHM can be still accounted for by orbital Doppler smearing 
if the mass ratio is larger than about 1/2. Again we consider the 
rapidly spinning white dwarf a less likely option as the white 
dwarf would be spinning at a velocity ~ 1/6 of the break-up ve- 
locity. 

Unfortunately we are unable to favor one or the other config- 
uration: case A and C have to be preferred if the phasing of each 
component is taken as "driving constraint"; while case B should 
be preferred if the measured FWHM (~ 1000 km/s) in the phase- 
resolved spectra is the most stringent requirement. In addition, 
none of the configuration is fully convincing as case A and C 
require a peculiar geometry of the hot-spot region, characterized 
by high thermal velocities and a "down-stream" temperature gra- 
dient. While case B "fails" in explaining the phase coincidence 
of the X-ray and ELR maxima. A fourth alternative configura- 
tion might locate the ELR and the optical pulse on the curved 
stream (accretion stream of flux pipe). This easily accounts for 
both the early appearance of the optical pulse and ELR with re- 
spect to the X-ray emission and for the phase coincidence of the 
X-ray and ELR maxima. However, none of the current models 
has yet explored such a possibility, hence the further discussion 
of the latter configuration is premature. 



At last, there are two more observational constraints imposed 
by the ELR. First the emission lines at ~4700 A is observed only 
in the phase-averaged spectra implying that it is very weak but 
visible across the whole period. Second, the line HeI/14859 is 
visible also at phase 0. 1 and at a higher significance than at phase 
0.3 and 0.5. These two facts are possibly evidence of a complex 
ELR geometry and/or of multiple line forming regions. 

5. Summary and conclusions 

The phase resolved spectra of HM Cnc presented in this paper 
have shown that its time resolved study is possible and in partic- 
ular it is feasible with FORS2 HIT-S mode off^ered at ESO-VLT, 
though care must be taken during the data reduction and the data 
analysis phase because of the low count level regime of the spec- 
tra and the deferred charges effect. 

Our phase resolved sequence has demonstrated the emission 
lines intensity varies across the period and that the emission lines 
in the phase-resolved spectra have smaller FWHM than those in 
the phase-averaged spectra. These are two strong observational 
constraints which however do not allow us -yet- to formulate 
any definitive conclusion about the location of the line form- 
ing region and on the nature of the period because of the low 
SNR and quality of our data. In particular, we cannot establish 
whether the emission lines arise from the X-ray emitting hot- 
spot or from the irradiated surface of the secondary star A higher 
SNR series of spectra would allow us to 1) establish whether the 
ELR is visible at any phase or they gets really "eclipsed"; 2) 
measure with smaller uncertainties the emission line centroids 
thus establishing whether they move across the period and the 
larger FWHM of the phase-averaged spectra is indeed an effect 
of Doppler broadening. Case A and C, which locate the ELR on 
the primary hot-spot should reveal emission lines with no o small 
Doppler shift. While, on the contrary, our case B configuration 
implies significant Doppler shift of the emission lines. All the 
discussed geometries are consistent both with unipolar induc- 
tor model and the direct impact accretor as those models do not 
prescribe any specific location of the ELR in their formulation. 
The face on-intermediate polar hypothesis, instead is not com- 
patible with the variable emissions Une intensity nor with their 
smaller FWHM in the phase-resolved spectra. Hence it does not 
fit any of the 3 cases discussed in section 4. In the double de- 
generate polar model, we expect that the ELR is on the primary 
and therefore it is consistent only with our configuration A or B. 
This model should be discarded in case large Doppler shift of 
the emission lines are measured. 

Finally, our phase-resolved spectral analysis showed that the 
Hell odd-term series emission lines at 4199. 8A and 4541. 6A 
were not detected (at any phase). This is consistent with the EWs 
inferred from the 2001 phase-averaged spectrum and confirms 
that the odd-term series lines are fainter than the even-term se- 
ries ones. The emission line centroids are not yet measured 
with sufficient accuracy to unambiguously address the issue of 
the presence of H in the spectrum of HM Cnc. However, we note 
that recent theoretical (D'Antona et al. 2006) and observational 
(Reinsch et al. 2007) studies of HM Cnc envisage the presence 
of H in DDB systems. 
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